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Abstract. We obtain convergent power series representations for Bloch waves in periodic high- 
contrast media. The material coefRcient in the inclusions can be positive or negative. The small 
expansion parameter is the ratio of period cell width to wavelength, and the coefBcient functions 
are solutions of the cell problems arising from formal asymptotic expansion. In the case of positive 
coefRcient, the dispersion relation has an infinite sequence of branches, each represented by a 
convergent even power series whose leading term is a branch of the dispersion relation for the 



p^ homogenized medium. In the negative case, there is a single branch. 

^ Key words: high contrast, double porosity, photonic crystal, negative index, dispersion relation, power 

'""5 series solution, generating function, Bloch wave, honiogenization, meta-material. 

in 

'TT ' 1 Introduction 



The objective of this work is the obtention of convergent power-series representations of Bloch waves and 
i -p, dispersion relations for the Helmholtz equation 



L^2 



V • ad(x)Vf7(x) + ^C/(x) = , (1.1) 



c 



, in which the scalar material coefficient a^i'x.) has small period d in each spatial direction and consists of two 

^ highly contrasting phases, with either positive or negative contrast 

^ ( 

•^ ±(Ph for x e d{P + m), m G Z", 

^ ad(x) = i 

fSJ I 1 for xe d(P= + m), meZ". 

r-^ The "inclusion" P is a subset of the unit cube Q in E" {n > 2) with C^ boundary dP, such that the 

^^ complement P'^ = Q\ P contains the boundary sides of Q (Fig. [I]). Both x and d carry units of length 

T-H so that x/d G K". The fixed constant 7 carries units of area, w is the operating frequency, and c is a 

►^ fixed reference celerity. We will develop power series in the expansion parameter d or the parameter of 

• 'H quasi-staticity, which measures the ratio of cell size to wavelength. 

/\ This equation has applications to acoustic waves in porous media, high-contrast photonic crystals, and 

S plasmonic crystals, in the regime of large wavelength-to-cell ratio. In the former, the inclusion contains 

the soft phase and the host material is of a stiff phase. In the latter two applications, the equation is a 
two-dimensional reduction of the Maxwell system in a periodic array of high-contrast rods, the field u is 
the out-of-plane component of the magnetic field, and ad(x) — e~^(x) is the reciprocal of the dielectric 
coefficient. The rods consist of a lossless dielectric material in the case of a positive coefficient or a lossless 
plasma in the case of a negative coefficient. 

As d tends toward zero in the positive case, a multi-branched homogenized, or quasi-static, dispersion 
relation relating frequency to Bloch wavevector emerges. This relation, which was obtained by Zhikov [HJ 
§8], [23^, involves the subset {/x„} of the set of all Dirichlet eigenvalues of the Laplacian —A = — V-V in P 
for which the means {4>n)p of the corresponding eigenfunctions 0„ do not vanish, 

fc^ = const4£^^I^. (1.2) 
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Figure 1: The unit cell in R^ with a high-contrast inclusion. 



Here, k is the wavenumber and the constant depends on the direction of the wavevector. The relation is 
graphed in Fig. [2] It reveals a sequence of spectral bands [/^*,/^„] characterized by those points on the 
abscissa for which the right-hand side of (1.2 1 is nonnegative. 



In the two-dimensional application to electromagnetic fields, the /x„ correspond to resonances of an 
effective magnetic permeability of the homogenized medium, as shown by Bouchitte and Felbacq [2 [S]. 
Modes that are localized in defects within a high-contrast medium of this kind have been shown to converge 
spectrally by Kamotski and Smyshlyaev |11| as well as Cherdantsev [1] using homogenization techniques. 
Limits of spectra in high-contrast media have also been obtained when the value of a in a periodic thin 
grating is small [71 [6] or when a is small in the periodically dispersed regions enclosed by the grating [22] . 

The remaining Dirichlet eigenvalues {/x^} of —A in P, that is, those whose eigenfunctions all have mean 
zero, also play a role in the small-d limit of equation ( 1.1 ). Hempel and Lienau [TU] proved that the spectra 
of the "double porosity" operators Ad = —V • ad(x)VC/(x) converge, as d — >■ 0, to the union of the bands 
[^* , fin] plus those eigenvalues /x^ that fall within the gaps between these bands. (In [TU], the authors actually 
deal with a fixed period and a material coefficient a that is fixed in the inclusions P+m and tends to infinity 
in the host material; by a scaling of a and x, the spectrum is seen to be equal to that of A^;.) Zhikov deals 
directly with the double porosity operator and identifies the limit of the spectra as ti —>■ in the sense of 
Hausdorff with the spectrum of a limiting operator A in the sense of two-scale convergence [251 §3] , where 
the two scales are x and y = x/d. 

At values oi ( = ^oj^/c^ belonging to a branch of relation (1.2), the equation {A — ()u — admits 
plane-wave solutions in x with a computable periodic microscopic variation in y e Q which is constant and 
nonzero in the host P'^ and variable in the inclusion P. On the other hand, at any eigenvalue ( = ^e or 
C = Mfi (^ ~ C)''^ = hs-s solutions that are supported in the inclusions alone and vanish in the host. If 
fi'^ falls in a gap of the homogenized dispersion relation, one views it as a band that has degenerated to 
a point, and if /i^ falls in a band of the homogenized dispersion relation, one views it as a gap that has 
degenerated to a point (see |10] for a discussion of this). Generically, eigenfunctions have nonzero mean, but 
for the important example of a circle (as illustrated in Fig. [2]) or a sphere, all but the radially symmetric 
eigenfunctions have zero mean. 

The guiding principle in the positive case is this: 



For each wavevector kk (\k\ ~ 1), there exists a unique number (^ ~ ^uP' j c? G \ii*^^ii„ 



the 



-,th 



spectral 



band defined by (1.2). If (^ is not a Dirichlet eigenvalue for —A in P, then, for rj = kd sufficiently small, 
the structure supports a Bloch wave with wavevector kk and square frequency uP' = c^CJ^ /^. As rj ^ Q, 
C,^ converges to (, and the Bloch wave converges weakly to a plane wave with wavevector kk and square 
frequency uP — c^C,/^ and has a computable strong two-scale limit. The Bloch wave and the frequency 
CP at nonzero values of rj are given by convergent power series in rj. On the other hand, for 
each Dirichlet eigenvalue ji^ for —A in P for which all eigenfunctions have mean zero, and for each kk, the 
structure supports a Bloch wave whose (scaled) square frequency C is given to leading order by ji'^ as rj — )■ 0. 

The demonstration of the convergent power series representation of the Bloch fields and the frequency, 
for both positive and negative coefficient in P, is the subject of this work. The case of negative material 
coefficient in P is simpler, as the homogenized dispersion relation has only one branch (Fig. w| and the 




Figure 2: The dispersion relation of the homogenized medium for an array of circles of radius 0.375 with positive 
material coefficient (t^ = 'jk'^ and (^ = ^ijp" jc? are the reduced square wavenumber and frequency). The asymptotes 
occur at the eigenvalues ^„ of —A that correspond to eigenfunctions Jo{finr), the spectral bands are the intervals 
[lj.^,fi„] where the function is positive, and the dots indicate the other eigenvalues n'^. 



eigenvalues of P do not play a role. We present all details for the case of positive coefficient for which the 
leading term Co of the power series for C'' hes in a band of the homogenized dispersion relation and does not 
coincide with /i^. In the case that Co = f^'e for some i, we allow the remark at the end of section p^ to suffice. 

Our approach to this problem is to expand the frequency and field in r] and then to prove convergence 
of the resulting power series to a dispersion relation for the structure and the associated Bloch fields for 
nonzero values of rj within some radius of convergence. The analysis extends our recent work on high- 
contrast plasmonic inclusions [HI |9] , in which we establish convergent power series representations for fields 
and dispersion relations in sub-wavelength plasmonic crystals (dispersions of rods consisting of a plasma 
material). An advantage of the method is that it does not rely on coercivity of the two-scale operator, and 
one is therefore able to prove existence of solutions for structures with large negative contrast in the material 
coefRcient. The advantage of a power series solution is that it provides an approximation of the true solution 
up to any algebraic order in 77. For certain coercive problems, it has been shown that the formal power- 
series expansion is an asymptotic series for the solution jTHl [2] and even that suitable truncation yields 
exponentially small error |12| . 

Convergent series for fields in composites with highly contrasting complex conductivity have been ob- 
tained by Bruno [31. As far back as 1890, Hermann Schwarz [16 obtained a power series solution of the 
equation Au+Xpu = in a bounded domain, where p{x) > 0. The expansion parameter is A, and the solution 
has a radius of convergence equal to the first eigenvalue. In fact, he also represents the first eigenfunction 
at the radius of convergence by a convergent sequence. 

Our strategy is as follows. 

1. Fix a wavevector kk and expand the field and frequency in power series in rj — kd to obtain an infinite 
coupled sequence of PDEs for the coefficients. 

2. Solve the first few equations to obtain the homogenized (quasi-static) dispersion relation and fields. 

3. Choose a branch of the homogenized dispersion relation, which corresponds to a unique value of C, 
and iteratively solve the infinite sequence of PDEs for the coefficients of the power series. In the 
positive case, it can observed explicitly how solubility of the system is obstructed by secular terms in 
the expansion if C is a Dirichlet eigenvalue of —A in the inclusion P. 

4. Obtain a nonzero radius of convergence of the series by establishing an exponential bound on the 
coefficients. This is achieved through the use of systems of generating functions. The radius is uniform 



over k in the negative case, and it is uniform in the positive case if C is bounded away from the 
eigenvalues ^^. 

5. Prove that the power series solve the field equation and that the odd coefficients of the series for the 
dispersion relation vanish. 

A remark on proof of convergence and comparison with plasmonic crystals. In [^ , we obtained power- 
series solutions for the dispersion relation and the associated Bloch waves for two-dimensional plasmonic 
crystals, in which the plasma frequency Wp tends to infinity as the inverse of the cell width. Following 
the scaling Wp — c/d of Shvets and Urzhumov [18], the small rj regime leads to a large negative dielectric 
coefficient in the cylindrical plasmonic inclusions ep = 1 — 1/(^^77^), where ^ — uj^ /{c^k'^). In that work, 
we use fine properties of the Catalan numbers to obtain quantitative values for a lower bound on the radius 
of convergence [21 §6bc] by obtaining an explicit exponential bound on the norms of the coefficients of the 
expansions. Special properties of the two-scale limit of the Bloch fields that were available in the plasmonic 
case but not in the case of positive contrast also played a crucial role in the analysis. The method of generating 
functions that we use in the present work is simpler and applicable as well to the case of plasmonic crystals, 
although it does not provide a way to estimate the radius of convergence. This approach can in fact be 
applied quite generally when one is able to obtain recursive bounds on the solutions of the higher-order cell 
problems, as, for example, for arrays of micro-resonators |17j . In [S], we assumed symmetry of the inclusion 
under a rotation of 180°. We do not make any symmetry assumption in the present paper, and our analysis 
here eliminates the need for this assumption for plasmonic crystals as well. 

2 Formal power series and dispersion relations 

We begin by fixing a wavenumber and introducing power series expansions of a Bloch field and its frequency 
in the parameter 77 = dk. Solving for the first two coefficient fields and then imposing the solvability condition 
at the next order in the matrix material leads to the homogenized dispersion relation. 

2.1 Long-wavelength Bloch waves 



The Helmholtz equation (1.1) can be written as the system 



flpc AxC/ + ^C/ = in the matrix, 

flpAxf/ + '^U = in the inclusion, (2.3) 

apc Vxf/|pc • n = apVyJJ\p ■ n on the interface. 

We investigate Bloch- wave solutions — solutions that are periodic at the microscopic scale and modulated by 
a plane wave at the macroscopic scale, 

C/(x) = M(y)e*'''^-^, x = dy, 

in which u(y) has the unit cube Q = [0, 1]" as a period cell. Apart from the definitions of Op and Opc, the 
equations for u(y) explicitly depend only on the unit vector k, the ratio w^/(fc^c^), and the parameter 77, 

ape fc2(Ay + r]2ik ■ Vy - r]'^)u + Tj^ '^u ^ 0, y e P", 

apk'^{Ay + r)2ik-Vy-r]'^)u + ri'^^u^0, y e P, (2.4) 

apc(Vy + ir/k)u\pc ■ n = ap(Vy + irik)u\p • n, y G dP. 

At this point, one can consider different asymptotic regimes depending on how the coefficients are defined. 
In this work, we take ap<: to be unity and Op to be proportional to the area of the unit cell. 

Opc = 1, flp = — = —--, (2.5) 

7 'yk'' 



where 7 is a constant carrying the dimension of area. It is also convenient to introduce the non-dimensional 
wavenumber r and square frequency C by 



7fc2 



C = 7- 



which transforms the system (2.4) into 



r2(A + ?7 2iK-V - if)u + tfC,u = 0, 
(A + ry 2ik -V -rf)u^-C,u = 0, 



yeP, 



T^(V + i-qk)u\p<^ ■ n = ?7^(V + irjk)u\p • n, ye 9P. 



(2.6) 



2.2 Power series expansions 

By inserting the power series ansatz 



u'' = Mo + lui + rfu2 + 
C" = Co + ryCi + ^^2 + ■ ■ 



(2.7) 
(2.8) 



into the system (2.6), one obtains a system of coupled partial differential equations for the coefficients, 



t2(Au„i + 2i«;-Vu„i_i -M„_2) + I]fco 0">n-2-£ = in P^, 

Au,„ + 2ik-VUm-l - Um-2 + Z^fcO Ci'Uji 

T'^(Vu„i + ikum-i)\p'= ■ n = (Vum_2 + iku 







m-3j|p 



in P, 
on 9P, 



(2.9) 



in which m,„ = for to < 0. For the function uq in the matrix P'^, (2.9) yields the BVP 



Auo = in P^ 

Vuq • n = on 9P, 



from which we infer that uq is a constant, which we denote by ug, in P"^. The system ( 2.9 ) yields the following 
boundary-value problem for uq in P: 



Amo + Couo = in P, 
uo I p = uo on dP. 



(2.10) 



At this stage (2.101 presents us with an alternative: either Co coincides with a Dirichlet eigenvalue, in which 



we take mo = if Co = M^ for some £; otherwise uq 7^ 0. In the former case, whether <^o feUs within a spectral 
band or a spectral gap, one obtains a power series for a Bloch wave for each wavevector (see the remark 
in section 3]) for which the leading-order term in the expansion of C'' is given by /i^. In the latter case, 
the possible frequencies Co &re determined by kk. At this point, we make the assumption that Co is not an 
eigenvalue and that uq ^ 0, which means that we are seeking a power series for a field that nonzero in the 
matrix. 

It is convenient to work with the dimensionless fields '(/',„ defined through 

Um = i™woV'm, 



so that -00 = 1 in P*^ and the system (2.9) becomes 



t2(AV'„ + 2K-VV'm-l + V'm-2) - Z^fcO {~iYQi'm^2-i =0 in P", 
(A + Co)V'm + 2k-Vlpni-i + V'm-2 + Z]"ll(~*)^C^ i'm-i = in P, 



T^(VV'm +KV',„-l)|i 



-(VV'm-2 + KV'm-3)|i 



on i9P. 



(2.11) 



Notice that the equation for ipm in P is a Dirichlet boundary- value problem for the operator A + Co with 
boundary data equal to the trace of ipm in P'^ and interior forcing coming from tpi in P for i < m and Ci 
for £ < m. The equation for ijjm in P'^, on the other hand, is a Neumann boundary- value problem for the 
Laplacian with periodic conditions on the boundary of the cube Q and normal derivative specified by those 
of ij^rn-i in P'^ and -0^ in P for i < m — 2 and interior forcing coming from tjjf in P'^ for £ <m and Q for 
£< m-2. 



The system (2.11) yields the following boundary- value problem for ^q in P 

AV-o + Co^o = in P, 



(2.12) 
tpo\p = 1 on dP. 



(2.13) 



This problem has a unique solution if C^q is not a Dirichlet eigenvalue of —A in P. The equation for ipi in 
P' is 

AV'i =0 in P", 

S/tpi-n + k-n — on dP, 
which has a unique solution subject to the mean-zero condition 

The equation for ■02 in P'^ is 

t2(A02 + 2k-VV'i -f l)-Co = in P^ 
T^(V-02 + 't0i)|p'^-n = -VVJolp-n on dP, 

which is a Neumann problem whose solvability is subject to the condition 

r^ / {k-Vipi + 1) - Co / V'o = . (solvability) (2.14) 

■/p<^ Jq 

This is implicitly the leading order of the dispersion relation, giving r^ as a function of Co (equivalently, k 
as a function of wq)- Indeed, the coefficient multiplying t^ is nonzero: By (2.13), we obtain 



J P" J P<: ^ J P" ' 



in which the inequality is strict because, since ipiiy) is periodic and k • y is not, Vi/'i is not a multiple of k. 
Thus, Jp, iVV'iP < |P1 and we obtain 

/ («;-V^i + 1) = / «-V0i + |P^| = - / |V0i|2 + |P-| > 0. 

J P" J P<: JP" 

Following Zhikov [21] . we can make the dispersion relation explicit by writing the solution tpQ in P in 
terms of the Dirichlet spectral data. 



oo 



My) = 1 + Co E ^^ My) = E ^""^ '^^.{y), y e P, (2.15) 

in which /z„ are those Dirichlet eigenvalues of —A in P whose i"-normalized eigenfunctions (/)„ appear in 
the expansion of the constant function in P and thus have nonzero mean. 



(0n)p = / 0«(y)rfy 7^0. 



The solvability condition (2.14) becomes 



« oo 



tJ-n{4>n) 



Co 



(2.16) 



As observed by Zhikov |21 
of values of Co such that 



this relation, shown in Fig.Mdefines an infinite sequence {Co ; '^ = 0, 1, 2, ... } 



*(") 



Mn < A^n < Co < A'n+l- 

The functions Q" ~ Co i'^^ ''') '-'^ '^ ^^'^ "^ define the leading order in rj of an infinite sequence of branches 
of the dispersion relation and reveals a sequence of stop and pass bands for a homogenized medium. A 
convergent series representation for C*-"-* for < 77 < Rnir), 

c(") = c^"^ + ^cf^ + ^'cr^ + ..., 

will be proved in section 4.2 subject to the condition that q" 7^ l^'i for all i. 



2.3 Negative material coefficient 

If we take the material coefficient in the inclusion to be negative, 

ap — 

7 






the plus sign before the Co in the system (2.12) becomes a minus sign, whereas the solvability condition 



(2.14) remains unaltered. The resulting dispersion relation is 



Jp- ^1 M« + Co 



(2.17) 



which differs from that for positive Up only by the plus sign in the denominator. The graph of this relation 
is obtained from Fig. [2] by refiection about the origin, and the result is that there is only a single dispersion 
relation that passes through the origin (Fig.^, as we require that ( = ^uj^ jf? be positive. 




Figure 3: The dispersion relation for the homogenized medium for an array of circles of radius 0.375 with negative 
material coefficient. AH positive values of C, are in the spectrum. 



3 Solution of the sequence of cell problems 

We shall develop the mathematics and state the theorems in the case of positive material coefficient in the 
inclusion, which is, in fact, the more complicated situation because of the resonances {/x^}. In section^ we 
discuss the modifications that need to be made in the case of negative material coefficient. 



When solving iteratively for the coefficients of the field and frequency, we do not assume a priori that 
the coefficients of the frequency are real. After proving that the series converges and solves the Helmholtz 
equation, we then prove a posteriori that the coefficients are real and that the odd ones vanish. 

The weak form of the cell problem for u is 



T / (Vw • Vf + irjk ■ (mVu — wVw) + 77 uv) — ?7 C / uv + 

J P-: JP'= 



+ rf / (Vii • Vw + irjk ■ {uVv - vVu) + rfuv) -rfC I uv ^ Q Vw G H^^^iQ), (3.18) 



and the weak form of the sequence of cell problems (2.11 ) for V'm is 



Jp" Jq 1,^0 

- / VV'm-2 -VV- [kV,„_3 • Vv - (k • VVm-3 + i^m-i) w] = V W G H^^^iQ), (3.19) 

in which ipm and Cm are set to zero for to < 0. 

Solving for ip^n in P. By restricting to test functions v with support in P and putting m i— 7> to + 2 in 



(3.19), one obtains a Dirichlet boundary- value problem for ipm in P, in terms of the boundary values of f/'n 



in P": 

/ (-VV'm • Vw + Co V'mW) - / [klprn-l • Vw - (k • VV'm-l + V'm-2) "] + 

, ■0m lap- = '0m|ap+- 
For TTi > 1, it is convenient to decompose V'm in P as follows (we define ipo = 0): 

■0m = V'm + (-O^CmV"*, 

in which i/im solves a Dirichlet problem involving ijjm in P'^, i't in P for i < m, and C« for £ < m, 

-VV'm ■ Vw + Co ■^m'i'j - / [kiJm-l • Vi) - (k • VV-m-l + V'm-2) «] + 

^ r?i— 1 
+ / V (-Z)'O ^m-^ t^ = V i- G ijl (^) , 

i'ni\ap~ = V'm|ap+- 
and f/;* is independent of to and satisfies 

(-VV'*-Vw + CoV'*i^)+ / V'oS = VveH^iP), 
p Jp 

V'*lap = 0. 
The strong form of this equation is 

(A + Co)0* = -V-o in P, 
0* = on dP, 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



and, given the expression (2.15) for "00, one can always solve explicitly for -0*, 



^* = J2 



\ (Mn - Co)' 



in P. 



(3.24) 



On the other hand, equation (3.22) for ipm can typically only be solved if Co 7^ ^^'g for all (., and this is 
precisely where this assumption plays its role. 



Solving tor ■;/;,„ in P^. The field %l;„t in P'^ is to be determined by (3.19) in terms of f/'f in P'^ for £ < m, 



ipe in P ioT £ < m — 2, and d ior £ < m — 2. Solvability of this equation is subject to two conditions. 



The first condition is obtained by setting u = in P'^, and is simply the problem (3.20) for '4'm-2 in P 



More precisely, given v G H}j^,j.(Q), decompose u = wq + wi according to the orthogonal decomposition of 

■ffper(Q): 

Hl,,iQ)=H^iP) + H^iP)^, 

where functions in Hq{P) are considered to reside in H^^^.{Q) by means of extension by zero. This means 
that wq = in P'^ and Aiii = vi in P. Now, given that V'm-2 satisfies (3.20) in P, we arrive at a Neumann 
problem for ipm in P"^ by replacing v with vi. 

The second solvability condition is obtained by setting v = vi = 1 in Q , 

{-^^-H.^-2 (J ^0 + Coy a) + I i](-*)'0^m-2-^ + 

+ Co / 0m-2 ~ T^ («;• VV-rn-l +^m-2) + / (k • VV-m-S + ^m-4) = 0. (3.25) 

Jp JP" Jp 

Given that this condition is satisfied, one can solve for ipm in P*^ up to an additive constant, and we will 
always take the unique solution that satisfies the zero-average condition 



1pm =0. 



(3.26) 



Solving for Cm- Equation (3.25) is to be viewed as the relation that determines the term Cm-2 of the 



power series expansion of the dispersion relation. Indeed, the quantity multiplying Cin-2 depends only on 
the inclusion, the values of r and k, and the choice of Co: ^nd it is always nonzero: 



f iPo + Cof iP* = \P'\ + f {1P0 + CoiP*) = \P1 + / £ T^^^^H^ </-« = \P1 + E 



Pl{(t>n) 



,tl (^» - ^")' 



> 0. 



'p Jp -^^i^ (^" ^ Co)^ 

Later on, we will prove that Cm = when m is odd. 

In order to solve inductively for all of the fields ipm and coefficients Cm , we will proceed from the established 
"base case" consisting of ?Ao in Q and Co , both of which have been specified above, as well as ipi and '02 in 



P"^. In P'^, -01 is determined by equation (3.19) for m = 1 



t2 / (V01 + k) . VV = VW e ifper(-P'), 

Jp= 



(3.27) 



IP 

i'l = 0, 

which has a unique zero-mean solution. Likewise, 02 is the unique zero-mean solution to 

T^ / (V02 + k)-Vv - T^ ( [k- V01 -I- l)v + Col ipov - I VtAo • Vt) == Vw e H^^^XP''), 
Jp" Jp'= Jq Jp 

^^2 = 0, 

(3.28) 

which is solvable because 0o satisfies (2.12) in P and because of the solvability condition (2.14), which the 
pair (r, Co) is assumed to satisfy. 



Theorem 1. //Co 7^ l^'t fof ciH i, there are unique functions ipm G H}^i,^{Q) and numbers C,„i G C such that 
''Pm = and Cm — for m < and such that 

1. the fields ipQ, ipi, and -02 *^ P'^ o,nd ipQ in P and the number Co coincide with those already defined; 



2. equations (3.191 and (3.251 are satisfied for all m; 
3- I "0™ — for all m; 



4- "0™ — 4'm + (—*)'" Cm "0* 'i'l^ Pj where i/),„ satisfies (3.22) 



Proof. The proof is by induction on the following statement that depends on ?i > 2: 
Statement Sn-' For n > 2, there exist 

• -0™ S Hp^^{Q) for TO < n — 2 with Jp^ -0™ ~Oifm>l, 
' 1pm & HleriP"") fom - 2 < m < u, with Jp^ 0^ = 0, 

• Ctrl G C for m < n — 2, 

that vanish for m < Q and coincide with ipo, 0i, and tp2 in P'^ , 0o in P , and Co as already defined, and such 



that (3.19) is satisfied for m <n and 0m — ipm + {—i)"^Cm'^* in P with 0„ satisfying (3.22). 

Statement ^2 holds by definition of "001 0i: a-nd ^02 in P'^, ipQ in P, and Co- Given the validity of Sm for 
m < n, we will prove Sn-^-l. 

Notice that (3.22 ), with to, = n — 1 is an equation for 0„_i with data that depend only on the functions 



"071 and numbers Cm already defined by the induction hypothesis. Since Co is not a Dirichlet eigenvalue of 
—A in P, this problem has a unique solution 'ipn-i- Let Cn-i be defined by equation (3.25) with to, = ri + 1. 
Besides Cm, that expression depends only on the numbers Cm and functions ipm presumed in the induction 
hypothesis. By definition of i/'*, the function defined by 0„_i := 0n-i + {—i)"'~^Cn-iip* satisfies (3.20) in 
P with TO = n — 1. Moreover, since ■0*|9p — and because of the boundary condition in (3.22), we have 
ipn-ildp- ~ V'n-ilap+i SO that the functions 0„_i in P'^ and P together define a function ipn-i that resides 
in H^^^iQ). 



Now consider (3.19 ) for m — n+1. The two conditions of solvability are satisfied. First, for all v e Hp^^{Q) 
with w = in P"^, this equation is satisfied by the definition of ipn-i- Second, by setting v — 1, one obtains 
the equation (3.25), which holds because of the definition of Cn-i- Thus (3.19) for m = n + 1 admits a unique 
solution tpn+i e H^^^iP") subject to Jp, V'«+i = 0. 



(3.19) 



To complete the proof, we must show that (3.25) holds for all to,. This follows from setting v = 1 in 

n 



Remark on interior eigenvalues 

In the case that Co — fJ'i for some £, we have seen that the higher-order cell problems in the foregoing analysis 
arc not solvable if it is assumed that mo ?^ in the matrix. Instead, one should set uq = in P"^, which is 
in accordance with the fact that the eigenfunctions for the two-scale limit vanish in the matrix and are a 
multiple of the eigenfunction for /z^ in the inclusion [22J . Now the leading-order term of the expansion of 
C is fixed at /i^, and one must solve for Cm ("t- > 0) and for all ■0m- Because Co is a Dirichlet eigenvalue 
in P, the solvability of ipm is subject to a Fredholm condition both in P"^ (where one sets u = 1) as well 
as in P (where one sets v = (/)). The solution is given up to an additive constant Pm in P'^ and up to the 
addition of 7m0 in P. The numbers /3m, 7m, and Cm are determined by subsequent solvability conditions. 
The procedure for solving for these constants and the fields '0m inductively turns out to be quite complex, 
though not insurmountable, and we believe that an exponential bound can be obtained using the technique 
of generating functions that wc illustrate in the next section. For this communication, we are satisfied to 
concentrate on the case Co 7^ A*^- 
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4 Convergence of the power series 

Convergence of the formal power series for u in H^{Q) and ^ in C is equivalent to an exponential bound on 
the norms iJ^-norms of the functions "ipm and the numbers Cm- An exponential bound is proved in Theorem^ 
below. The iterative scheme for solving for ipm and Cm provides recursive upper bounds on their norms. 
These bounds are based on the standard H^ estimates obtained in the appendix. 

4.1 Bounds on the fields 



From (3.19) and the bound (7.62) in Problem 2 of the appendix, with K = i7pc max{l. A\, we obtain 



r'^Wi^TnWHHP-) < K 



{^Wlllm-lWrniP") + \\i^Ta-2\\m(P-')) + \\An-2\\m(P) + 



+ 2||V'm-3lli/i(^) + IIV'm-4|ki(P) + ^ IC^ (II V'm-2-£ ||i/i (P=) + || V'm-2-£ || //i (P) ) 



1=0 



(4.29) 



From equation (3.22) for -0 and the bound (7.60), we obtain 



IV'mllHi(P) < Kr 



||V'm|lHi(P<^) +2||V'm-l|lHi(P) + |lV'm-2|lifi(P) + ^ ICd II V'm-fH H1(P) 



1=1 



(4.30) 



in which K^ < Cit^ + C2 if Co is bounded away from the eigenvalues /i^. 

In (3.25), the constant multiplying Cm-2 is at least IP'^I for all r, and we obtain, by adjusting K if 
necessary. 



ICml < K 



CollV'mllHi(P) +T^ (|lV'm+l|ki(pc) + llV'mllffi(P-)) + 

771—1 

+ ||'0m-l|Ui(P) + ||V'm-2||_H-i(P) + ^ IC^I {\\i>ra^i\\m{P) + llV'm-^? || ffi (P<^) ) 



t=l 



(4.31) 



From the decomposition ifj^ — "ipm + (~i)™CmV'* and the bound ( |7.58 ) on -0*, we obtain 

IIV^mllHi(P) < IIV^mllHi(P) + S^ICml, (4.32) 

in which Br = (Sit^ + Sa)^. 

4.2 Proof of convergence 

The foregoing analysis is most transparent using the quasistaticity 77 as the expansion parameter, as we have 
done. For the investigation of structures of fixed cell size with varying wavenumber k, it is more suitable to 
expand the fields and dispersion relation in powers of the cell size. To this end, set 

d rj 



Then we have the expansions 



Recall the relation Um = i"^uoipri 



h T 



Uq + prui + p^t'^U2 + 

2„2, 



C = Co + prCi + pV"C2 
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Theorem 2. Let a branch Co = CoC^j''") G [/^nuA'm) of the homogenized dispersion relation he given. For 
each r > such that Co 7^ M^ V^, t/iere exist rea/ positive numbers C and J , independent of k, such that 



\\r"'^m\\HHQ)<C.r 



\T'^Cm\<cr 



(4.33) 



For each tq > and for each e > 0, C and J can be chosen such that (4.33) holds for all k and for all r 
such that < r < To and 

Co e [K„, Aim) \ \J{p-'i - e, /4 + e). (4.34) 



If the band [/i*^./im) contains none of the /.i'^, then C and J are independent of e. 
Proof. Let us simplify our notation for the frequency and the norms of the fields, 

Pm = |k"Vm||Hi(P<=), 

Pm = |k"Vm||ffi(P), 

Pm = ||-r"Vm||Hi(P), 

'^TTi I ' Sm I ■ 



The bounds ( |4.29|4.30|4.31|4.32| yield the following bounds for m > 1: 



(rn — 1 

(rn— 1 
Pm + 2tP,„_i +T^p„_2 + Y. SiPm-i,^ 
e=i ' (4.35) 

?n— 1 \ 

TPm+l + T^p^ + Tp„j_i + t'^P„i-2 + L sKPrn-*" + Pni-f) ) > 

^=1 ^ 



S^ < K 



[CoPn 



B s 



To prove that the numbers Pm, Pm, Sm, and Pm are exponentially bounded, it suffices to prove that the 
sequences of positive numbers, defined by the system of recursion relations below, obtained by replacing 
inequality with equality in the relations above, are exponentially bounded. Indeed, one observes that pm < 
Om, Pm ^ bm, Sm < Cm, and Pm < ^m • Recall that we have defined ipo — 0, which gives po = 0. 



Om + l = K[2Tdm + T^flm-l + dm-l + 2r(i„ 



Omi Ojyi Cjyi Ojjyi U 

flo =Po, ai =Pi, 6o =Po = 0, Co = |Co|, ^o =Po, 

7n — 1 
-2+T d„i_3 + ^ Cg[dm-l-l + dm-l- 

J2 cedm-e ) 
i=i ' 

(m—l 
Co bm + TClm+l + r'^am + rdm-l + T^dm-2 + J2 Ct{dm-t + dm-t) 
i=l 

O'm Oyji 



bm = Kr[dm+ '^rdm-l + T^dm-2 



JD-yCj) 



m < 0), 

m> 1), 

771 > 1), 

771 > 1), 

771 > 1). 



(4.36) 

In addition, it is proved in the append ix tha t B^ = [Bit^ + ^2)^ (equation 7.58) and that, under the 

condition (4.34), Kr < Cit^+C2 (equation 7.61). This means that, if ri < T2, the values of {pmiPm, Sm,Pm} 



for T = Ti are bounded above by the values of {a™, bm,Cm, dm} for r = T2 (observe that the initial conditions 
are independent of r). To prove the theorem, it therefore suffices to prove the sequences {0^, bm, Cm, dm} 
are exponentially bounded for each r > 0. The bound can be taken to be independent of k by choosing 
0-1 > Pi for all k. 

It is convenient to transform the recursion relations above by the change of variables 
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The system for {om, bm, Cm, dm} is equivalent to the following system for {am, bm, Cm, dm}'- 

O-m — bm+l ~ Cm+1 = rfjn+l = 

flo ==_po, ai =Pi, bi^po^ 0, ci == ICol, di ^ po, 

/ m— 1 

a-m ^ K(2Tam-l + T^am-2 + dm-1 + '2Tdm-2 + T^dm-3 + J2 Ci{am-l-i + dm-t) 
^ 1=1 

bm = Kr ( O-m^l + 2r(i„i_i + t'^ dm-2 + Z) Qf^m+l-f) ) 
^ 1=2 ' 

( m-1 \ 

Cm = -fi^l Co fem + •ram + T^a„j_i + rd„_i + T'^dm~2 + E Q(am-f + drn+l-f) 

V «_9 / 



m < 0), 

TO > 2), 

TO > 2), 

TO> 2), 
TO > 2). 



(4.37) 



We shall prove that {«„, 6^, Cm, d,„} are exponentially bounded by proving that their generating functions 
have a nonzero radius of convergence. Set 

oo 

Y, dnz'' = za{z) + ao, 

n=0 

CO 

E &„^" = zP{z), 

CXD 

X; c„2" = Z7(z), 

n=Q 

oo 

X; dn^" = z5{z), 



(4.38) 



The system (4.37) is equivalent to the following system among the formal power series q;(z), fi{z), 7(2), and 
5{z) 

z{a{z) - ai) = K{2Tz'^a{z) + t22:^(zq;(z) + oq) + (z + 2x2^ + t2z^)z5(z) + 

+ z{z^{z)){za{z) +ao) + (z7(z))(z(5(z))j , 
z(/3(z) - 61) = X,(z2a(z) + (2rz + t''z'')z6{z) + 2-1(2(7(2) - ci)){z{5{z) - di))) , 
z(7(z) - ci) = /i:('Tz(a(z) - fli) + Co z{l3{z) - ai) + T^z'^a{z) + {tz + t^z^)z5{z) + 

+ (z{-i{z) - c^))za(z) + z-Hz{^{z) - ci))(z(<5(z) - di))) , 
z(,5(z) - di) = z(/3(z) - fei) + S,z(7(z) - ci). 

Define the following functions of five variables: 

A{a, 13,^,5, z) = -(a-ai) -+ K{2Tza + T'^z{za + qq) + {z + 2tz^ +t^z^)6 + 

+ zj{za + ao) + z"fS] , 
B{a, 13, 7, (5, z) = -(/3 - 61) + K^ (za + (2rz + t^z^)S + (7 - ci)(<5 - di)) , 
C(a, /3, 7, 5, z) = -(7 - ci) + K(T{a - ai) + Co(/3 - ai) + r^za + {tz + t^z^)5 + 

+ 2(7 - ci)a + (7 - ci)((5 - di)j , 
Z?(a, /3, 7, (5, z) = -(<5 - di) + (^ - &i) + B,(7 - ci). 

One can check that A, B, C, and Z? all vanish at (a, /3, 7, 5, z) = (ai, 61, Ci, di, 0). Moreover, the determinant 
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of the Jacobian matrix with respect to the first four variables is nonzero at this point: 



d{A,B,C,D) 



d{a,l3,-f,6) 



z=0 



-10 

-1 KriS-di) Kr{l-Ci) 

Kt KCo -l + K{S-di) K{j-ci) 

1 Br -1 



and thus we obtain 



cLet ^l^'^J^'? (ai,bi,ci,di,0) = 1^0. 



(4.39) 



By the implicit function theorem of analytic functions of several variables, we infer that there exist analytic 
functions represented by power series 

a{z) — ai + a2Z + a-iZ^ + . . . , 
/3(z) = 6i +622 + ^32:^ + •••, 

7(2) = Cl + C2Z + CgZ^ + . . . , 

5(z) = di + diz + d'iZ^ + . . . , 
that are convergent in a neighborhood of z = 0, say for |z| < i? > 0, and such that 

(AS,C,7^)(a(z),/3(z),7(z),J(z),z) = (0,0,0,0) 



for \z\ < R. By definition of A, B, C, and D, the functions a, /?, 7, and S satisfy ( |4.38 ), and their coefficients 
therefore satisfy (4.371. By the convergence of the power series, the coefficients and therefore also Om, 6m, 
Cm, and dm are exponentially bounded. D 



5 Power series solution of the cell problem 

We are now ready to prove that the formal power series in fact converge to Bloch wave solutions and their 
dispersion relations for small but nonzero r/, as long as (^Q{k,T) is not a Dirichlet eigenvalue of A in P. In 
the process, we also show that the odd coefficients Cm vanish. 

5.1 Solutions for nonzero rj 

The exponential bound for the functions T™-ipm and the numbers r™Cm and therefore also for the functions 
Tpm and numbers Cm, implies that the formal power series for the field u and the square frequency C converge 
for \r]\ < R for some R > 0. We show that the functions defined by these power series are indeed solutions 
of the cell problem, for each r > 0. For |?7| < R, define the functions 

u^ — uq + rjui + ifu2 + . . . , 

C" = Co + Ci'? + C277' + • ■ • , 

in Q, in which Um — i™uoV'mj and V'm are the solutions to the cell problems already described. For 
V G Hp^j.{Q) and rj < R, define 



a^iv) 



(Vu'' • Vv + irjk ■ {vPWv - vWvP) + jfvPv) - jfC / W^v + 



+ Tf j (Vu'' • Vw + irjk ■ {vPVv - vVvP) + ifu'^v) - rfC, / u''u. (5.40) 
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This function of 77 has a convergent power series in rj that is obtained by inserting the series for u^ and (^ 
into the expression for a''(u) and expanding in powers of j]. The coefficients of this expansion are exactly the 
right-hand side of equation ( |3.19 l, muhiphed by z™, with the identification Um = i™MoV'm- But this is equal 
to zero for all m because the functions i/'m satisfy (3.19). Thus all coefficients of the power-series expansion 
of a''(w) vanish. We conclude that a^{v) — for all v G Hp^j.{Q) and jr/j < R, which means that (u^X^) 
satisfies 



which is the weak formulation (3.181 of the PDE for u, for all \ri\ < R 



Theorem 3. In the solution {vP , C), all of the functions ipm, where Um — i^UQipm, cire real-valued, Q„i G 
for all m, and C^m = if ni is odd. 



Proof. From the weak form of the cell problem (3.18), one shows that C, is necessarily real-valued by putting 
V = u^ and therefore all coefficients Cm are real. One then proceeds by induction on the following set of 
statements for n > 2: 

1. Tprn in P is real- valued if to < n — 2, 

2. 1pm in P"^ is real- valued if m < n, 

3. C.m — ^i "ni < n ~ 2 and m is odd. 



These statements hold for n = 2, as one observes from equations ( 2.15|2.13|2.16|3.28 ) that ipo in Q, ipi and 
■02 in P'^, and Co are all real-valued. Let n > 2 be arbitrary, and let us prove these statements for n replaced 
by n -|- 1. By equation (3.22), with to = n — 1, we see that V'n-i is real-valued in P. Equation (3.251 with 



n -I- 1 shows that (- 



Cn-i is real (?/'* is real-valued) . Since Cn-i is real, C»i-i = if n — 1 is odd. 



Thus ipn-i — V'n-i + (— 0" ^Cn-iV'* is rcal-valucd. Finally, equation ( |3.19[ ) for to = n + 1 shows that ipn+i 
is real- valued in P'^. D 



5.2 Convergence of as 77 —> 

The solution of the cell problem yields a Bloch wave solution for our periodic structure. 



[/"(x) 



CXJ 



Mn+l 



(kx/rj) 



which tends, in the sense of strong two-scale convergence, to a separable function that is a plane wave in x: 

U''^ix)^Uo{x,y)^e'''^--Uoiy). 

In the two-scale theory [2T],[22j §3], one writes the solution as a function of x plus a function of {x,y) that 
is supported in the soft phase P, 



Uo{x, y) = un{x) + i;o(x, y) = e 






+ e 



Co 2^ 



f-n — Co 



(Pniv), 



in which the sum over n is extended by zero into P^ and we keep in mind that Co depends on k and k. 
The function ^0(2;) is a Bloch solution in a homogenized medium. This homogeneous medium is realized 
approximately by a real periodic structure with small cell width d, in which the fields are measured in 
the matrix only, and the role of the inclusions is to lend the bulk medium its rich band-gap structure. 
This idea is discussed in several other works such as [U [H [HI HOI [HI HI] j in which the extreme properties of 
periodically dispersed inclusions hosted by a neutral matrix material produce meta-materials with interesting 
bulk properties. 
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6 Negative material coefficient 



If we take the material coefficient inside the inclusion P to be large and negative, while retaining a small 
positive index in the matrix, 

apc — I, Op — — T-^, 

7 7fc^ 

one can follow the minus-sign through all of the foregoing calculations, obtaining a convergent power series 
for the field u and frequency w, with analogs of all of the theorems we have proved. The situation is actually 



simpler because now the equation in P is always solvable and the nonresonance condition (4.34) is no longer 



needed. We show in the appendix (see equations 7.53 7.57 7.59 1 that, in place of the constants K, K^, and 
Br in section |4.2[ we may use a single constant K that is independent of k and r. We therefore obtain a 
radius of convergence that is valid for all wavevectors. 



In place of equations (2.12) and (2.15) for the field inside P, we obtain the Helmholtz equation with a 



change in sign in front of the frequency Co: 



A^o - Co^o = in P, 
V'oU = 1 on dP, 



and the solution is 



My) - E 



Mn {4'n)p 



0«(y), y e P. 



(6.41) 



(6.42) 



The expression (2.14) for the solvability for ^^2 in P'^ remains unaltered, and, using the explicit expression 
for ipQ, yields 

^J■7^{(|)n)p 



t' f («;.V^i + l)-CoE 



\ M« + Co 



(6.43) 



Since the field i/"! in P'^ is decoupled from the inclusion, it does not depend on apc^ and therefore this 
dispersion relation differs from that for positive ap<: only by the plus-sign in front of Co in the denominator. 
It possesses therefore only one branch, which passes through the origin in (r, Co)-space, or equivalently, in 
(fc, a;)-space, as illustrated in Fig. [3J 

The weak form (3.19) of the sequence of PDEs in the unit cell becomes, for negative index, 






\7v + T^ [ktpjn-1 • Vw - (aC • VV^m-i -I- 1pm-2) «] + / E (~*)'^Cf ■0m-2-<? V + 

Jp^ Jq ^^0 

+ / V^„_2 • Vi; + / [kV,„_3 • Vw - (k • VV'm-s + V'm^4) v] = VvE H^^^iQ), (6.44) 

and the equations for ip, ip, and -0* in the decomposition ipm — ipm + (— *)™Cm'0* in P become 

/ (VV'm • Vu -I- Co V'rnW) + / [k?A„i_1 • Vw - (k • V'f/'m-l + V'm-2) v] + 

^ m 

+ / y^i-iYCiM^iv - Vve Hl{P), 

^mlap- = V'mloP+J 

^ VMi • Vw -I- Co -^mWJ -I- / [klp„i-l • Vu - (k • VV'm-l + i>m~2) v] + 

^ rn — 1 

Jp »_1 



(6.45) 



.46) 



i=i 



i^nildP- = MldP+i 
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and 



(Vip, • Vw + Co Av) + tpov ^ V w e i?o (P), 
p Jp 



(6.47) 



^*\dP = 0. 
The strong form of the latter equation is 

(-A + Co)V'* = -i^o in P, 
V'* = on dP. 

and the solution in terms of the Dirichlet eigenfunctions of P is 

fJ-n {4>n)p , 



^*--E 



in p. 



- (A^n + Co)2 "" 

The solvability condition for f/'m in P^, which determines the value of Cm- 2, is 



(6.48) 



{-l) Cm-2 



i'O + Co MJ*)+ V (~l)^C<? V'm-2-f + 

VJQ Jp / JQ ^^;^ 

+ Co / '0m-2 ^ "^^ / ('«• V?A,„-1 + V'm-2) - / (k • V-0m-3 + '/'m-4) = 0. (6.49) 

Again, the quantity in multiplying Cm-2 is always nonzero, 



hpo + cJ A ^ |P1 + E 
Jo Jp 



^i(a*„ + Co)^ 



> 0, 



so that one can always solve for Cm.-2- 

7 Appendix 

In this appendix, we derive the bounds on the fields in P and P'^ that are employed in the recursive system 



of inequalities in section 4.2 



7.1 Bound for the resolvent of the Dirichlet- A in P 

Let {ipj, Vj^'^i, < I'l < z/2 < i^3 < • ■ • , be the Dirichlet spectral data for the region P, where the Lpj form 
an orthonormal Hilbert-space basis for L'^{P), that is, 

—Atpj = h'jfj in P, 

(fij =0 on dP, 

W'^jWl^p) = 1- 

The set {vj} is the disjoint union of the sets {/in} and {/i^}, where the former consists of those eigenvalues 
for which there is an eigenfunction with nonzero mean. Using integration by parts, one finds that V(/)j is 
orthogonal to V(/)j' in the mean square inner product if j ^ j' and that 



Now let u be the solution of 



W^'PjWl^p) = ^j\\'Pj\\h(p) = ^j- 

Au + vu^G in P, 
u — Q on dP, 



(7.50) 
(7.51) 
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in which v is not an eigenvalue and G E L^{P). We shall obtain a bound on the solution u in H^{P) in 
terms of G. 

If the eigenfunction expansion of G is 



G^^bjipj, 



then 



j=o 



— - — (jC,, and Vu = >^ — - — Vw,-, 
j=0 ■> j=0 ■> 



so that using (7.50) we obtain 






and |lVu|l|2(p) =^i^j 

3=0 



Hence, 



l"llL2(P) + l|Vu||i2(p) 









< max 



l + Uj 



] \\v- Vj 






Since the latter sum is equal to ||G||?2/d\, we obtain 



L2(p), 



in which 



||'"||h1(P) < C',y|lG||i2(p), 

r (l + '^.O^ 

6,y = max ■' 



(7.52) 



3 W-l'jl 

Now let us bound the constants C±,^(, in terms of r. 

In the case of negative material coefficient, the pair (r, Co) lies on the homogenized dispersion relation 



(6.43), and we must put ly = — C^g in (7.51). Putting ly = —(q in (7.52) yields a uniform bound on C_^(,: 

C-Co < Co for Co > 0. (7.53) 



In the case of positive material coefficient, (t,Co) lies on the homogenized dispersion relation (2.16), 

■^ fJ'n{4'n)p 



BCo 






\ Mn - Co ' 



(7.54) 



in which B = [/pc(K-V'0i + 1)]^^. Let us choose the branch of this relation associated with the eigenvalue 
pLm for some fixed m. This means that Co(''") lies in the band [/i*^^,^,^) for all r, where /i^ is the root of 
the right-hand-side of ( 7.54[ ) that lies between ^im-i and /im for ttt, > 1 and /ip = 0. The right-hand side of 
(7.54) is positive for /i*„ < Co < ^^m and tends to infinity as Co tends to /i„i from below (see Fig. [2]). Let K 
be such that 



^Cc 



-^ tJ-n{(Pn)p 

J ^^n - Co 



<K, 



so that 



0<BC( 



IJ'm{(Pm}p ^ _2 I 7^ 



Mm — Co 



<t" 
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This yields the bound 



1 



Co 



< 



[SCoAi,n(0m>p] ' (r^ + K)< [B^i:^^lm{(bm)l] ' (t^ + K). 



(7.55) 



From (7.52) 



C'Co — T^T-OX ■ 



(1 + fJ'n 



Mr7 



Co 



,M 



in which 



^' + "'^'<M if^^M« 



li'i -Col 

Here is precisely where the issue of secular terms in the series solution arises. If there is an eigenvalue /i^ 
in [fj,'^, fj,^), then M — ^ oo as Co ~^ M^i ^^'^ ^^ ^ oo when C,q = fi'g. In this case, there is no uniform bound 
on C(^g for bounded values of t and, moreover, for that t for which Co = M^i the problem (7.51 1 is not well 
posed. 

In the case in which there are no values /i^ in [/i*„,/i„i), M can be taken to be a fixed number for all 
Co G [^J'm^ Mm) and thus for all r. More generally, if small neighborhoods of those eigenvalues /x^ that are in 
IfJ,'^, iJ-m) are excised from this band, M can be taken to be finite. Thus, using our estimate for Cq^ together 
with (7.551, we obtain the bound 



Cco < Ait'^ + A2, Vr such that Cq(t) G [^;;,, /i^) \ |J(/x^ - e,^^ + e). 



(7.56) 



in which Ai and A2 depend on P and the branch of the dispersion relation, and, in case there is an eigenvalue 
/i£ in [/i^, /im)i also on e, and A2 — > 00 as £ — > 0. 

7.2 Bounds for the fields in P and P^ 



The field xp^. In the negative case, we obtain, from the explicit expression (6.48) for ip, 



1^* 



Lffi(P) 



E 



(1 + ^i^)^il (0„)^ 



(Mn + Co)^ 



<s, 



(7.57) 



which is bounded by the single number B for all Co ^ 0. 



In the positive case, the expression (3.24) of V'* yields 

(l+/i„)/l2 ((/.„)p 



WMh 



HP) 



E 



(^" - Co)4 
Since all of these summands except for n = m are uniformly bounded over n and Co £ [Mm7/^ni)j ^^ obtain 

(1+ Mm) 5 Mm (0m)p 



IIV'*IIhi(^) < k + 



{^■■m - Co)^ 



in which K depends only on P and the branch of the dispersion relation. In view of (7.55), we obtain, for 
some constants Bi and i?2. 



Mhmp) <{B^T' + B2Y=:Br. 



(7.58) 



Problem 1: Field extension from P^ to Q and solution in P. Suppose that ip is defined as an H^ function 
in P'^, G e L^{P), and C is not a Dirichlet eigenvalue of —A in P. Consider the extension of jp to f{^{Q) 
such that 



I {\/ij-\Iv~Qi!v)^ I Gv VveH^iP), 
ilj\dP+. 



IP 

MdP 
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Decompose ?/; as "0 = V-'o + "01 such that 



/ (V01 • Vw + Vi«) = VveH^iP), 



and 



We obtain 



IP 

i^ildp- ='4'\ap 



(V0O • Vz; - CV-o^) = / (G + (1 + OV-i)" yv € iJoH^) 
V'olap- =0 



p Jp 



II^iIIhi(p) < const • llV'll^i(gp) < MHmiP-), 
UoWh^p) < CciWGWmP) + \l + CIUiWhHp)), 

in which, as we have discussed, C^ is large when ( is close to a Dirichlet eigenvalue of —A in P. From these 
two bounds, we obtain 

UWrniP) < Gcl|G|U2(p) + AiCc\l + CI + l)\mHHP^)- 
Now, in the case of negative material coefficient, when (r, Co) lies on the homogenized dispersion relation, 



we must put C = —Co- Because of (7.53) we obtain a bound 



Mhhp) < KiWGh^fp) + UWH^Pn) (7-59) 

that holds for all k and t. 

In the positive case, we must put C = Co, where (r, Co) lies on a chosen branch of the dispersion relation, 
and we obtain 

Mhhp) < Kr{\\G\\mP) + ||V^||,,i(pe)), (7.60) 

in which 

K^ = maxjGco, A(G(;Jl + Mml + !)}■ 



If the branch contains none of the fif or these numbers are excised, the bound (7.56) yields 



Kr<CiT^+C2 Vr such that Co (t) e [//;;, Mm) \|J(M£-£, A^^+e)- (7-61) 

e 

in which Ci and C2 again depend on P and the branch, and, in case there is an eigenvalue /i^ in [/ij„,/im), 
also on e, and C2 — >■ 00 as e — >■ 0. 

Problem 2: Solution in P'^. Suppose that the function 1}^ E Hp^j.{P'^) satisfies 

/" V?A • Vw = / (Fi • Vv + Gi«) + f {F2-S/V + G2V) Wve Hl^^iQ), 
Jp^ Jp^ Jp 

Tp = 0. 



Ipa 

The following solvability conditions are necessarily satisfied. 



f iF2-Vv + G2v) = WveH^iP), 



Gi + / G2 = 0. 
P" Jp 
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The zero-average condition implies the Poincarc bound 

MmiP") < ^p-||VV'||l2(pc). 
Extend tp to Hp^j.{Q) such that Af/; = ■i/' in -P so that 

and set v — ip above. 

||VV'||i2(pe) < (lli^illL^(pc) + ||Gi|U2(p.))||i/.||Hi(p.) + {\\F2\\l-(p) + \\G2\\mp)) \\i>\\m{P) 

< (||Fi|U2(p.) + IIGilU.(pc) + A (WF^h^P) + ||G2J|l^(p))) Mh^p^)- 

Using the Poincare bound, wc then obtain 

Mhhp') < ^Ic (||i^i||L2(pc) + ||Gi|U2(p.) + A(||F2|U2(P) + ||G2||l2(p))) . (7.62) 
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